WSU-HEP-0805 

MCTP-08-44 

NSF-KITP-0861 



c ■ 

1—5 ! Gordon L. Kane^, Alexey A. Petrov^'^, Jing Shao^ and Lian-Tao Wang^ 

^ I ^Michigan Center for Theoretical Physics, 

I"!' University of Michigan, Ann Arbor, Michigan 48109 

"^Department of Physics and Astronomy, 

Oh; 

QJ ', Wayne State University, Detroit, Michigan 48201 

I ^Physics Department, Princeton University, Princeton, New Jersey 08544 

(N ■ 

> ■ (Dated: June 11, 2008) 

: 

^ ; Abstract 

'""i ' In principle particle spins can be measured from their production cross sections once their mass is 

in ■ 

. approximately known. The method works in practice because spins are quantized and cross sections 
^p. ' depend strongly on spins. It can be used to determine, for example, the spin of the top quark. Direct 

■ application of this method to supersymmetric theories will have to overcome the challenge of measuring 
mass at the LHC, which could require high statistics. In this article, we propose a method of measuring 
the spins of the colored superpatners by combining rate information for several channels and a set of 
kinematical variables, without directly measuring their masses. We argue that such a method could 
lead to an early determination of the spin of gluino and squarks. This method can be applied to the 
measurement of spin of other new physics particles and more general scenarios. 



PACS numbers: 



I. INTRODUCTION 



The upcoming experiments at the Large Hadron Colhder (LHC) will shed light on the solutions 
of the shortcomings of the Standard Model (SM), and should identify new degrees of freedom 
associated with the New Physics. Most of the leading candidates for such physics involve super- 
symmetry, which predicts doubling of the observed particle spectrum. It is expected that some of 
the supersymmetric partners of the SM particles will be produced at the LHC. It will be essential 
to confirm that the candidate superpartners indeed differ by half a unit in spin from their partners 
in order to be confident it is indeed supersymmetry that is being discovered. A rigorous determi- 
nation of the spins of new particles will be difficult and complicated, particularly early with limited 
statistics, and will require a detailed analysis of decay angular distributions and correlations 
Though it will be necessary to do this eventually, we want to argue that an initial determination 
of the spins can in many cases be made, particularly for particles that dominate production. If 
the new physics is indeed supersymmetric the initial production of colored superpartners is likely 
(though not certain) to dominate, and we focus on a preliminary measurement of their spins. As 
an illustration of the method, we apply it to the determination of the spin of the top-quark, whose 
measurement by this approach does not seem to have been previously reported. The top-quark 
analysis illustrates both the strength and limitations of the method, clearly distinguishing a spin-| 
"top quark" from a spin-0 one, but not from a spin-1 one. 

A standard way of testing a theory, such as supersymmetry, is by comparing theoretical pre- 
dictions with experimental data. Given preliminary observation of a signal for a new particle, the 
minimal test would simply involve calculation of the production rate for a particle of a given spin, 
color, and other quantum numbers, and checking that the theory and data are consistent. For 
instance, the production rates of colored superpartners, like that of the top quarks, are largely de- 
termined by their QCD interactions, so these checks do not depend significantly on other features 
of the chosen model of physics beyond the SM. 

Such a test can only be performed once the particle's mass is measured, since the production 
rate depends on the mass. In order to test a particular model, one can in addition compare the 
production rates computed within this model with that of some alternative scenario - comparisons 
are normally less sensitive to systematic uncertainties such as cross section normalization errors. 
For example, for the top quark, one can compute the rate for the particle in the same representation 



of color group, but of different spin, say zero or one. These rates would differ significantly both 
because of the number of spin degrees of freedom and because of different angular momentum 
configuration of the final state; depending on their spin, the final state particles can be produced 
in either s— or p— wave. As it turns out, the latter effect is significant even at a hadron collider. As 
we shall see, even an approximate measurement of the mass provides good discrimination between 
the cases of different spin, simply because the resulting production rates are very different. This 
implies that one can determine the spin of the produced particles even with the limited-luminosity 
early measurements of production cross sections and mass. The results can of course be improved 
with more data and analysis. Similarly, we will compare alternative spin scenarios for color octet 
"gluinos" . The method can later be extended to test both the color and spin separately. 

Of course accurate calculations of cross sections of strongly-interacting particles require inclusion 
of higher order QCD effects. In order to demonstrate the effectiveness of the method, here we only 
consider the rates computed at the tree level. Indeed, the choice of the renormalization and 
factorization scales, as well as other effects, will bring additional uncertainty into evaluation of the 
absolute rates. For a given mass, the effects of such uncertainties will be essentially the same for the 
different spin cases. However, as we shall show in this paper, the relevant comparison is between 
different spin scenarios with different mass scales. Therefore, the effect of higher order QCD 
corrections do not cancel in this comparison. Nevertheless, we shall argue that such corrections 
will not change the rate in ways that significantly affect the effectiveness of our method. 

It can be challenging to measure the mass of a particle, particularly at a hadron collider. 
Most importantly, typical kinematical variables, such as Hj- or rries, essentially measure the mass 
differences between the produced particle and its decay chain members. Absolute mass scale can 
seldom be determined from these variables in scenarios such as supersymmetry, where the end 
product of a decay chain is a massive neutral stable particle which escapes detection. With only 
information about rates in certain channels, and mass differences, as we will discuss in detail below, 
there are typically degeneracies between supersymmetry and other scenarios with different spin 
partners. 

Studies of mass determination using detailed kinematical information were performed in 
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211 ]. Application of these methods typically requires rather 



high statistics. Moreover, such mass measurement would have to be fairly accurate in order to 
be useful for spin determination. Since the rate scales with a high power of mass, a small error 



in mass measurement could translate into a large uncertainty in rate predictions. In this paper, 
we suggest a rather different method. Instead of trying to measure the mass directly, we propose 
to combine the rate of several channels and a set of commonly used kinematical distributions. 
We demonstrate that this method allows us to break the degeneracy and could lead to an early 
determination of the spins of the new particles produced at the LHC 

Ideally, one would also like to measure gauge couplings, color charges, and electroweak charges 
of the new particles, which are also predicted by the theory, and also test spins of other particles 
(e.g. charginos and neutralinos) that are produced with smaller rates and/or in decay chains. 
Such rates and branching ratios do depend on particle spins, so once the spins and masses of the 
particles that are dominantly produced are determined, some tests or determinations of others 
may be possible, depending on how well the relevant signals can be isolated. Some constraints 
on color charges can be checked fairly simply. For example, if gluinos are produced from color 
octet gluons, their color representation can be8®8 = 1 + 8 + 8 + 10 + 10 + 27, and one can 
calculate the production cross sections holding other quantum numbers fixed. It is appropriate to 
proceed by initially checking one property at a time (such as spin) , setting others (such as color) at 
their expected values. As more statistics become available, several variables can be simultaneously 
varied. 

This paper is organized as follows: we motivate our discussion by considering measurement of 
the top-quark spin in Sec. II, provide basic set-up in Sec. Ill and apply the proposed method to 
determination of the spin of gluino in Sec. IV. We conclude in Sec. V. 

II. MOTIVATION: DETERMINING THE SPIN OF THE TOP-QUARK 

In this section, we briefly discuss the method of determining the spin of a candidate top-quark, 
denoted "top-quark," with production rate information. The production of the top quark at the 
Tevatron is dominated by the s— channel process qq tt, which implies that the cross section is 
proportional to the velocity (3 of the produced top-quark near the kinematical threshold. However 
if instead scalar particles are produced, the configuration of the final state must be CP even and 
so the cross section is dominantly p— wave and is proportional to near the threshold. Thus we 
expect the spin-^ "top quark" to have a larger production rate than spin-0. In Fig{Tl we plot the 
cross sections as functions of mass of the "top quark", where we also included the gluon fusion 
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FIG. 1: Cross section as a function of the "top quark" mass at the Tevatron. The lower dashed hne is 
for the case of spin-^ "top quark", and the sohd hne is for the case of spin-0 "top quarks". The two 
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process. The production cross-section for spin-0 "top quarks" is indeed significantly below that for 
the spin- 1 "top quark". Thus given the mass of the "top quark", the magnitude of the production 
-section strongly favors the case of spin-| "top quark" ^. 
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Ill. CROSS SECTION FOR COLORED PARTNERS WITH DIFFERENT SPIN 



In this section, we compare the production cross sections of colored superpartners and their 
counterparts (with different spin) in alternative scenarios. 

For the purpose of illustrating the idea, we implemented a simple same-spin scenario^ where 
for each SM quark q there is a massive partner q'. Results for models with more extended fermion 



^ Indeed, it should be remarked that if a top-quark candidate decays into a SM b-quark and a W boson, it can 

only have a half-integral spin. 
^ One should note that this is not the case of the Universal Extra-Dimension (UED) scenario 27 1. 



partner sectors could be obtained by scaling from our result. The coupling relevant for production 
is schematically, 

l^int = q'L^'qL + u'ji^'uR + d'ji^'dn (1) 

We begin with spin (gs)^, \ ( q) and 1 color octets {g'). The QCD interactions of these 
color octets with the quarks and gluon are completely fixed by gauge invariance. T he p arton cross 



sections for pair production of these particles are already calculated in reference 
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or exam ple, at the LHC, gluon-gluon fusion leads to the following elementary cross-sections 
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These results are obtained including diagrams with s— channel gluon fusion and t— channel gluino 
exchange. For the spin-0 and spin-1 cases, there is also a four-point interaction diagram. From 
the above result, we see that they have the same functional dependence on (3, so the argument 
applied for the top-quark case does not apply here. However the differences in the cross sections 
originate from the differences in the interaction form determined by the spin structure. As for 
the top quark case, one again finds the longitudinal enhancement for the cross section of the 
spin-1 gluino. In Figj2l cross sections at LHC of the three cases are plotted as functions of the 
mass. In our calculation, the factorization scale fip and renormalization scale fin are set to be 
fJ-F = fJ-R = '^g,g',gs- From the plot, we can see that the cross section increases by roughly an 
order of magnitude as spin increases from to 1 for a given mass. On the other hand, these cross 
sections decrease rapidly with mass. Roughly speaking, for a increase in mass by 200-300 GeV, 
the cross section decreases by an order of magnitude. Thus to match a given result of the cross 
section, the masses of gluino for different spin must be different by 200-300 GeV. This makes it 
possible to determine the spin by a rough mass measurement. 

For later reference, we also present here the numerical results for squark pair (fermionic partner 



^ We consider a real scalar here. 




FIG. 2: A plot of the cross sections for gluino pair production (solid line), spin-1 gluon partner g' pair 
production (dashed line) and spin-0 gluon partner gs pair production (dot-dashed line) at LHC. In the 
calculation, extra color triplets (e.g. q or q') are taken to be 5 TeV. 

q' pair) production and squark-gluino {q'-g' ) associated production in Fig. [3l For the associated 
production, the cross section depends on both the squark {q') mass and the gluino {g') mass. To 
give a simple example and for later convenience, we choose gluino mass and g' mass such that 
their pair production rate are matched. Then we plot the cross section as a functions of squark 
(g') mass. 



The cross sections are calculated using MadGraph/MadEvent 40|]. In the calculation, we use 
the parton distribution function CTEQ6L [41]. The renormalization scale and the factorization 
scale lip are chosen fJ.R = IJ-f = and IJ.r = fip = 0.5 x {rrigg/ + rriggi) for squark (or q') pair 
production and squark (g') gluino {g') associated production, respectively. 

Our calculation of the parton cross sections is done at tree-level. The next-to-leading order 



calculation usually changes the result by a factor of order unity 30, ImI, l32] • As the multiplicative 
K-factors always increase the production cross section, such corrections should not change the 
overall hierarchy of the rates for different spin assignments. The enhancement factors of the cross 
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FIG. 3: Top: Cross sections for squark pair production (solid line) and fermionic quark partner q' pair 
production (dashed line) at LHC. In the calculation, gluino and g' are taken to be 5 TeV. Bottom: Cross 
sections for squark-gluino associated production (solid line) and q'- g' associated production (dashed line) 
at LHC. The gluino mass and g' mass are 608 GeV and 850 GeV respectively in order to match the cross 
sections of gluino pair production and g' pair production. 

section for particles with different spin but in the same color representation have the same color 
factor and the same initial state. Although not identical, those factors are not expected to be 
significantly different. Notice also we are not comparing cross sections at the same mass scale. As 



we will discuss in detail below, we will match the cross sections of particles with different spin with 
different masses. However, those different masses are typically within a factor of 2. Therefore, 
additional effects from different choices of renormalization and factorization scale will have to be 



included. We note that such differences in K-factors are typically quite small, about 10% [301. l32| 
across the mass range we are interested in. We have also partially taken this into account by 
chosing /iR and /x^? to be correlated with the masses of the particle produced. 

IV. SPIN DETERMINATION 

In this section, we present our study of spin determination for the possible cases where de- 
generacy could occur. First, we briefly review and clarify the degeneracy by combining rate and 
mass differences. For concreteness, we will focus on the comparison between supersymmetry and 
the same spin scenario defined in Eq. [H We will assume the final decay products always include 
stable neutral particles. Of course, this would be the LSP in the case of SUSY. We will assume 
the existence of such a stable particle, denoted by A, in the case of same spin scenario due to the 
implementation of certain discrete symmetry. Examples which have qualitatively the same feature 
as such a scenario include Universal Extra-Dimensions {27] and little Higgs models with T-parity 
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For simplicity, we will only consider simple decay chains following the production. A study of 
a set of models with more complicated decay chains is beyond the scope of this initial paper. We 
will briefly comment on such scenarios in the conclusion. 

A. Degeneracy with rate and mass differences 



As pointed out in Ref. 38|] and [6|], typical transverse variables are largely sensitive to the 
mass differences between the initial and subsequent particles in the decay chain. This is also a 
physics reason for the set of degeneracies in measuring mass parameters within supersymmetry. 



pointed out in Ref. [39|. Therefore, in a particular channel, there can be a degeneracy between 
supersymmetry and an alternative scenario with the same-spin partners. Schematically, suppose 
the production cross section for some superpartner of mass m is cr^^^^. We could choose the mass 
of a similar partner in the alternative scenario, m', so that the production cross section matches. 



^mi^ ~ '^'rn^^ ■ Same time, we could learn from some transverse kinematical variables 

the information about mass differences in supersymmetry, say Am = m — mLSP- However, if 
those variables mainly contain information about mass differences (which is true to a pretty good 
approximation), we can adjust the mass difference in the alternative scenario. Am' = m' — tua, to 
fix the kinematical distribution as well. If other information is available about the mass, that can 
remove the degeneracy, but typically that will not occur with early collider data. 

We demonstrate this degeneracy in a simple example. Consider the comparison between a 
squark and some fermionic quark partner q'. Let's also assume they are both pair produced, and 
followed by similar simple decay chains q —>■ q + LSP and q' q + A, where A is some stable 
neutral particle, analogous to the lightest supersymmetric particle (LSP) in SUSY models. In 



our simulation, we generate parton-level events with MadGraph/MadEvent [40|] for the squark 
pair production and decay pp —>■ qq* ^ qq + 2 LSP, as well as the q' pair production and decay 
pp q'q'* — > qqAA Then the parton-level events are passed to PYTHIA 6.4 [42] for parton 
shower, fragmentation, and PGS4 [43] for detector reconstruction. It is straightforward to match 
the production rate by choosing appropriate rrig and nig'. For example, m^ = 549 GeV and 
rriq/ = 900 GeV both give rise to a ^ 2 pb. In the above examples, we have not decoupled the 
color octet particles, but this will not affect our discussion below. 

Next, we adjust the mass of the missing particle to match the kinematical distributions. Notice 
that choosing the mass difference to be the same. Am = m^ — mLSP = i^q' — f^iA) will give 
rise to significant difference in kinematical distribution, particularly when the mass difference is 
large. Indeed, in order to match the rate, we have to choose rriqi to be significantly greater than 
rriq. Therefore, choosing a large and fixed mass difference usually implies Am ^ ^t^lsp- In fhis 
case, the available energy ~ Am tends to be shared evenly between the LSP and the rest of the 
observable decay products {q in this case). On the other hand, in the alternative scenario, we 
usually have Am much closer to the missing particle mass m^. It is therefore typical for the 
recoiling quark to have a harder spectrum. In the FigJH we plot the Ht ^ distribution for q and q' 
pair productions. In both cases, we keep the same mass difference Am = 452 GeV, i.e., mLSP = 97 



Here, for simplicity, we haven't include the qq and q' q' final states, and only include the quark partners of ur in 
the final state. 
^ Effective mass is defined as Ht = fir + ^ Pt 
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FIG. 4: Left: distributions for squark pair production (solid line) and q' pair production (dashed 
line) at LHC with both the cross section and mass difference matched. The mass parameters involved 
are ruq = 549 GeV, m,/ = 900 GeV, mLSP = 97 GeV and uia = 448 GeV. Right: Ht distributions for 
the same processes but with the mass difference adjusted in the same spin scenario. The mass of the 
missing particle now is = 548 GeV. In both cases, we have generated 3000 events, corresponding to 
approximately 1.5 fb~^ integrated luminosity. 



GeV and = 448 GeV. One can check that the difference in the average of the Ht does match 
correctly with the theoretical expectation. 

However, this does not imply the absence of the degeneracy. The key requirement for the 
existence of such a degeneracy is that the transverse variables are approximately functions of only 
mass differences, even if the functions may be different for different scenarios. In this case, it is still 
possible to adjust the mass difference in both scenarios to obtain indistinguishable distributions, 
even though the mass differences will be different in those cases. Although this requires very 
specific masses, it is important that we address this possibility in order to have a robust spin 
determination method. We can see this again in the example of squark production and decay. For 
the examples we just considered, if Ma is increased to 548 GeV, then one cannot distinguish the 
Ht of these two cases as seen in the FigjH 

A similar conclusion applies to three body decays such as ^ ^ gg + LSP. More explicitly, 
we consider supersymmetric and same-spin partners, which have about the same production rate 
(cr ~ 0.8 pb) and mass splitting. 
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FIG. 5: Left: Ht distributions for gluino pair production (solid line) and g' pair production (dashed line) 
at LHC with both the cross section and the mass difference matched. The mass parameters involved are 
m-g = 800 GeV, m^/ = 1060 GeV, rriLSP = 137 GeV and tua = 397 GeV. Right: Ht distributions for 
the same processes but with the mass difference adjusted in the same spin scenario. The mass of the 
missing particle is now = 497 GeV. In both cases, we have generated 3000 events, corresponding to 
approximately 3.75 fb^^ integrated luminosity. 

• mg = 800 GeV, ttilsp = 137 GeV 

• rUg, = 1060 GeV, = 397 GeV. 

The Ht distributions for both productions are plotted in Fig. O Again we can adjust the mass 
difference in the same spin scenario such that the two distributions match. This can be seen in 
FiglSl where we have increased niA to 497 GeV. 

Before proceeding further, we remark that there is a possibility that we could tell that the part- 
ners are from the same spin scenario and not supersymmetry. Assume the new physics discovered 
at the LHC is actually in the same spin scenario. We want to verify whether a supersymmetric 
scenario could match the rate and the same kinematical distributions. Because of the difference in 
production rate, the supersymmetric scenario will necessarily have a much lighter colored partner, 
squark or gluino, or both. Then, in order to macth the Ht distribution, we must have at least 
comparable mass gaps between the colored partner and the missing particle. Therefore, there is a 
possiblity in this case that we cannot find a solution for the supersymmetric scenario. For example, 
from the previous discussion, we know that, to match a cross section about 0.8 pb, either a gluino 



with mass rrig = 800 GeV or a g' with mass rUg' = 1060 GeV is required. However, if the mass 
difference is measured to be larger than 700 GeV, then the supersymmetric scenario is excluded 
since otherwise the LSP mass is below the LEP II bound. 

B. Matching mult i- channels 

As we have seen so far, there are purely kinematical differences for models with different mass 
of the decaying particle. We also notice that these differences are generally not very significant for 
measuring spins since the mass difference Am can be adjusted independently of the mass of the 
decaying particle. A more detailed analysis for these situation is certainly necessary to explore such 
differences for spin determinantion. However, in reality, there will be more handles to distinguish 
spin using the rate information. A complete model usually contains several other new particles 
beside the color octet In many theories, we expect the masses of the colored partners to be similar, 
i.e., around TeV scale. Therefore, at least several of them will be copiously produced at the LHC. 
We should consider the production channels of those colored states together. We will then have 
more observables which we could use to break the degeneracy discussed above. In this section, 
we present a method to achieve this goal by combining information from several channels. We 
show that this method, based on a set of simple observables, is effective and independent of mass 
measurement. 

To demonstrate our method, we focus on distinguishing supersymmetry and the same spin 
scenario. We start with the case in which the squark mass is lighter than the gluino mass; so 
gluinos decay through a two-body process. In our example, the gluino and squark masses are 
taken to be nig = 608 GeV and rriq = 549 GeV. Then we find that to match the gluino pair 
production rate, the g' mass is fixed to be mg> = 850 GeV. Furthermore we match the gluino- 
squark production rate by varying the q' mass, which leads to m^/ = 780 GeV. This can be seen 
from Fig. [6l where the two vertical dashed lines correspond to the masses of g and q' while the 
horizontal dashed line corresponds to 20 pb for the cross section of the associated production. 
After this matching, the cross sections of qq and q'q' productions are fixed to be about 8 pb and 
19 pb respectively. The cross sections for both models are calculated in MadGraph/MadEvent and 
are 

SUSY: agg ~ 4 pb, agq ~ 20 pb, a^g ~ 8 pb 




FIG. 6: Production cross sections for gq production (solid line), g'q' pair production (dashed line), qq 
production (dash-dot line) and q'q' production (dash-dot-dot line). The gg and g'g' production cross 
sections are already matched by choosing gluino mass rrig = 608 GeV and g' mass mgi = 850 GeV. The 
horizontal dashed line corresponds to 20 pb for the cross section of the associated production. If we match 
to this rate, the squark mass is 549 GeV and the q' mass is 780 GeV as indicated by the vertical dashed 
lines. After this matching, the cross section of qq and q'q' productions are fixed to be about 8 pb and 
19 pb respectively. 

Same Spin Scenario: ag'g' ~ 4 pb, cr^/g/ ~ 20 pb, cr^/g/ ~ 19 pb (5) 

One can see that the rate of q' pair production cannot be matched with that of the squark pair 
production at the same time as other rates are matched. 

More generally, we notice that the cross section of gluino {g') pair production is almost inde- 
pendent of the squark (g') mass. Therefore, by matching the cross sections for gluino and g' pair 
production, the g' mass is fixed for a given gluino mass. After that, the gluino-squark {g'-q') cross 
section only changes as squark {q') mass varies. Fig. [6] shows a plot displaying the cross sections as 
functions of squark {q') mass with gluino {g') mass fixed. As we can see from this figure, generally, 
the three production rates cannot be matched at the same time. In addition, such matching will 



typically force us to consider quite different mass splittings. Therefore, we also expect the typical 
jet transverse momentum variable, such as Ht-, will be quite different. 

The discussions above on matching individual cross sections are obviously a simplification to 
demonstrate the principles. In reality, after showering and forming the jets, the channels considered 
above can give sizable overlapping contributions to the same final state. Therefore, instead of 
matching the cross section for each channel, we perform next a more realistic study by matching 
directly to observables, in particular, the jet multiplicity and pt distribution ( again using Hq^ 
as a representative example). Generally we expect no degeneracy even by considering the jet 
multiplicities. The reason is as follows. Since we are including another particle in the analysis, 
there are three mass variables in the models. These three variables can be fixed by matching the 
total cross section and the two Ht peaks (if the masses of the two color particles are considerably 
different). Because of the difference in the spins in the two models, the ratios of the cross sections 
between different channels are generically different. Therefore the jet counts are not likely to be 
the same, because these channels have different jet-multiplicities. 

1. Case I: nig > nig 

Let us start with the case where the gluino is heavier than the squark in the SUSY model. For 
simplicity, we take it to be the SUSY model with mass parameters nig = 608 GeV, niq = 549 GeV 
and mLSP = 97 GeV. The total cross section of the gg, qq and gq final states is approximately 
31.6 pb. Then by varying the mass parameters nigi and nig' in the same spin scenario, one can 
match the total rate of the SUSY model; several choices of masses are listed in Table HI For each 
of these models, we generate 2000 events in MadEvent-BRIDGE-Pythia-PGS setup^ and compare 
the jet multiplicities with those of the SUSY model. The decay of the gluino is a two-body process 
g ^ q + q, followed by the decay of squark q q + LSP. Here we include the first two generations 
of squarks for simplicity^. For the same spin scenario, we consider the similar decay processes of 
g' and q': g' ^ q' + q and q' ^ q + A. Once the jet counts are obtained, the difference between 



^ We use BRIDGE[4J| to decay the parton events containing the g\\imo{g') and squark((7') before Pythia and PGS. 
^ The third generation squarks((7') will decay into b-jets which can be tagged and studied seperately. 



SUSY and same spin scenario can be characterized by defining a x-square-like quantity 

1 N /^SUSY _ SamcSpin \ 2 

In the following, we consider jet counts for = 5 bins with rijet = 1,2,3,4,5+. Here nf^^^ 
^^SamcSpm^ is the number of events in each bin for the SUSY (same spin) model. The standard 
deviation is defined by erf = (^o-SameSpm<j2 _|_ ^^susy^2^ where o",, = ^/rii + 1 for each model. In our 
simulation, we have used the PGS4 default detector configuration and cone jet algorithm with a 
cone size 0.5. In addition, we use the PGS trigger with low thresholds 

• Inclusive fir 90 GeV 

• Inclusive single-jet 400 GeV 

• Jet plus (180 GeV, 80 GeV) 

• Accoplanar jet and (100 GeV, 80 GeV, 1 < < 2) 

• Accoplanar dijets (200 GeV, Acf) < 2) 

To be realistic, we also impose the selection cuts for jets: Pt > 50 GeV and r] < 2.5. 

In Table m the average of Ht and the value of AS* are shown for given nig', uiqi and tua- For the 
SUSY model in the comparison, the average of Ht is 913 GeV, and the jet counts are {189, 965, 
620, 171, 46} for jet number Uj = 1,2,3,4,5+. As can be seen from the table, values of AS" are 
generically large for cases with m^/ > m^/ . This can be understood since in this case q' and g' decay 
very differently from squark and gluino. However, for cases with m^/ < nig/, AS can be smaller. 
The same spin model with the smallest AS* in the Table is the one with rUg' = 900 GeV, niq' = 800 
GeV and ttia = 550 GeV, where AS = 2.5. The jet multiplicity and Ht distribution are shown in 
FiglTl First of all, one can see that AS is not small, and so the jet multiplicity is not completely 
matched.^ Second, even if AS" is small, the Ht distribution is different enough to eliminates 
the possible degeneracy. So based on the above observations, one can see that the degeneracy is 
unlikely to exist considering total cross section, Ht and jet multiplicity simultaneously. 

^ It actually cannot be further reduced by tuning the mass parameters. 



ruA (GeV) 


Same Spin Models (m^/, niqi) (GeV) 


(1000, 640) 


(950, 720) 


(900, 800) 


(850, 880) 


(800, 960) 


(750, 1120) 


100 


(1218, 6.6) 


(1271, 7.9) 


(1307, 7.2) 


(1373, 7.7) 


(1399, 13.4) 


(1474, 16.4) 


250 


(1104, 5.8) 


(1178, 6.4) 


(1218, 6.2) 


(1254, 8.0) 


(1289, 16.8) 


(1339, 20.0) 


400 


(875, 3.2) 


(984, 5.1) 


(1036, 4.8) 


(1057, 7.2) 


(1072, 17.2) 


(1103, 18.4) 


550 


(584, 6.5) 


(682, 2.7) 


(767, 2.5) 


(776, 3.9) 


(758, 14.6) 


(820, 12.0) 


700 




(452, 18.7) 


(427, 6.0) 


(410, 7.4) 


(456, 6.4) 


(706, 14.2) 



TABLE I: This table shows {Hx, AS") for each same spin model parameterized by three mass parameters 
nigi, ruq' and for Case I. Here AS is defined in Eq.([U]). In each column, a set of (m^/, m^/) is chosen to 
match the total cross section with that of the SUSY model, which is defined by rrig = 608 GeV, rriq = 549 
GeV and rriLSP = 97 GeV. In each row, a mass is specified, which can be seen from the first column. 
For the SUSY model in the comparison, Ht = 913 GeV, and the jet counts are given by {189, 965, 620, 
171, 46} for jet number rij = 1,2,3,4,5+. The data indicate that Ht and jet multiplicity of the same 
spin model cannot be matched to those in the SUSY model simultaneously 

2. Case II: rrig < nig 

For the case where the squark is heavier than the gluino in the SUSY model, we consider the 
following example: rrig = 640 GeV, rriq = 800 GeV and mLSP = 100 GeV. The total cross section 
of the gg, qq and gq final states is approximately 11.3 pb. In such a case, the jet multiplicity can 
be matched better for same spin models with q' heavier than g'. First we vary the masses of g' 
and q' to match the total cross section, where four possible choices of {mg',mq') are listed in Table 
HTF . The decay processes are simplified by taking q ^ g + q and g q + q + LSP, and similarly for 
q' and g'. As seen from the table, the best-fit of jet multiplicity is achieved when nigi = 950 GeV, 
niqi = 1120 GeV and = 600 GeV, where AS = 0.8. The jet multiplicity and Ht distribution 
are shown in FiglHl Though the jet multiplicities are well matched, the differences (both the 

^ The reason to choose these four cases is that m^' — m^' is close to rriq — rrig. For large mass difference, it should 
have larger kinematical difference. 
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FIG. 7: Left: The jet counts for the SUSY model (sohd hne) in the Case I and the same spin model 
(dashed line) which has the smallest AS* in Table I. We vary mass parameters rrig, rrig and ttilsp for the 
SUSY model, and rugi, niqi and niA for the same spin model. Both models are tuned to have the same 
total cross section. Right: The Ht distributions for the same pair of models: the SUSY model (solid line) 
and the same spin model (dashed line). Both plots contain 2000 events, corresponding to approximately 
70 pb~^ integrated luminosity. 



width and the position of the peak) in the Ht distributions should be enough to distinguish them. 
Therefore, again we see that matching the Ht distributions and matching the jet multiphcities are 
in conflict with each other and cannot be satisfied simultaneously. As long as squarks are not too 
heavy and the production rate is of about the same order of magnitude as that of gluino, the above 
considerations are reasonable. The analysis may be difficult when squarks become much heavier 
than gluinos. 

So far we have demonstrated that there is finally no degeneracy for the cases we studied. 
However, for realistic models, both the squarks (g') and gluinos [g') could decay through many 
different channels, leading to longer decay chain and more final-state particles. This could change 
the jet multiplicity and jet Pt significantly, and it is not clear whether there are degeneracies or 
not in these situations. If there were, one would expect including more observables would again 
eliminate the degeneracy. On the other hand, when LHC data are available, one may be able to 
learn something about the decay topology or decay branching ratios. Given that information, the 
potential model dependence can be reduced. Finally, if there is no degeneracy, then one can figure 



niA (GeV) 


Same Spin Models (m^/, niqi) (GeV) 


(1000,1020) 


(950, 1120) 


(900, 1240) 


(850,1400) 


100 


(1530, 9.6) 


(1657, 5.6) 


(1713, 7.2) 


(1692, 5.8) 


200 


(1478, 8.9) 


(1585, 5.5) 


(1643, 7.0) 


(1639, 6.1) 


400 


(1258, 10.3) 


(1358, 4.3) 


(1377, 4.6) 


(1381, 2.6) 


600 


(914, 13.6) 


(984, 0.8) 


(1015, 2.8) 


(1037, 5.6) 


800 


(474, 20.4) 


(560, 15.3) 


(701, 19.8) 


(961, 22.5) 



TABLE II: This table shows {Hx^ AS) for each same spin model parameterized by three mass parameters 
m^', rUqi and rriA for Case II. Here AS" is defined in Eq.([U]). In each column, a set of (m.^', rriqi) is chosen to 
match the total cross section with that of the SUSY model, which is defined by rrig = 640 GeV, rriq = 800 
GeV and ttilsp = 100 GeV. In each row, a mass ruA is specified, which can be seen from the first column. 
For the SUSY model in the comparison, Ht = 1161 GeV, and the jet counts are given by {7, 80, 387, 
689, 834} for jet number rij = 1,2,3,4,5+. The data indicate that Ht and jet multiplicity of the same 
spin model cannot be matched to those in the SUSY model simultaneously 

out the spin of gluino or squark by fitting to the data. 
V. CONCLUSIONS 

We present a new method of determining the spins of TeV scale new physics particles by 
utilizing production rate information. Basically, if the mass and the production cross section are 
measured, the spin is then determined. For many forms of new physics that might be discovered, the 
method will work well. Degeneracies may occur, in which case further analysis is needed, which is 
provided. After clarifying the degeneracy in spin measurement by using rate information in a single 
production channel, we propose to combine information, such as jet counts and Ht distributions, 
from several channels. This method will be particularly useful for an early determination of spin 
of the colored new physics particles since we expect at least several of them will be produced 
copiously. Such a method does not require precise mass measurement. In particular, we consider 
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FIG. 8: Left: The jet counts for the SUSY model (sohd hne) in the Case II and the same spin model 
(dashed line) which has the smallest A5 in Table II. We vary mass parameters nig, niq and mLSP for the 
SUSY model, and rugi, niqi and niA for the same spin model. Both models are tuned to have the same 
total cross section. Right: The Ht distributions for the same pair of models: the SUSY model (solid line) 
and the same spin model (dashed line). Both plots contain 2000 events, corresponding to approximately 
200 pb~^ integrated luminosity. 



two cases where the SUSY models are characteristically different: nig > rriq and nig < niq. For 
the first case with a specific set of SUSY mass parameters {nig, niq and mLSp), we find that the 
jet multiplicity cannot be well matched by the same spin model with similar mass parameters 
{nig', niq/ and m^). In addition, the distribution shows extra difference from that of the SUSY 
model. For the second case, we find that the jet multiplicity can be well matched by the same 
spin model. However, the difference in the Ht distribution can be used to distinguish these two 
scenarios. The reason why this method works relies on the fact that the production cross sections 
encode the spin information of the SM partners in a nontrivial way, which lead to signals which 
can be used to distinguish scenarios with different spin assignments. In some cases, the pure rate 
signal may be hard to see due to special mass splittings. However, the kinematical differences are 
generically present and allow us to distinguish these scenarios eventually. 

Detailed study of spin determination using our method is certainly needed. First of all, although 
we have imposed certain selection cuts which mimic those would be used in Standard Model back- 
ground suppression, a precise evaluation of the effectiveness of our method can only be obtained 



after carefully including full background and more realistic detector effects. 

We have demonstrated our method using simple decay chains, which could be the dominant 
ones in realistic models. At the same time, more complicated decay patterns could well occur. On 
the one hand, such decay channels, especially those with leptons, will greatly enhance the discovery 
potential and signal to background ratio, which could improve the sensitivity of of method. On 
the other hand, they could introduce more model dependence. A detailed study of disentangling 
those effect is certainly important and useful, which will be discussed in follow-up papers. 
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